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Abstract. Data from e^e~ annihilation into hadrons, collected by the JADE experiment at centre-of-mass 
energies between 14 GeV and 44 GeV, are used to study moments of event shape distributions. Models 
with hadronisation parameters tuned to the LEP 1 precision data provide an adequate description of the 
low energy data studied here. The NLO QCD calculations, however, show systematic deficiencies for some 
of the moments. The strong coupling measured from the moments which are reasonably described by NLO 
QCD, 

as(Afzo) = 0.1287 ± 0.0007(stat.) ± 0.0011(exp.) ± 0.0022(had.) ± 0.0075(theo.) , 
is consistent with the world average. 

PACS. 12.38.Bx Perturbative calculations - 12.38.Qk Experimental tests 
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1 Introduction 

Electron-positron annihilation into hadrons constitutes a 
precise testing ground of Quantum Chromodynamics 
(QCD). Commonly jet production rates or distributions 
of event shape variables have been studied. Predictions 
of perturbative QCD combined with hadronisation cor- 
rections derived from models have been found to describe 
the data at low and high energies well, see e.g. [2-6]. 

In this analysis we use data from the JADE experi- 
ment, recorded in the years 1979 to 1986 at the PETRA 
e+e~ collider at DESY at six centre-of-mass (cm.) ener- 
gies y/s covering the range 14-44 GeV. We measure the 
first five moments of event shape variables for the first 
time in this low ^/s region of e^e~ annihilation and com- 
pare the data to predictions by Monte Carlo (MC) mod- 
els and by perturbative QCD. Moments sample all phase 
space, but are more sensitive to specific parts of phase 
space, dependent on their order. From the comparison of 
the data with theory we extract the strong coupling as. 
The measurement of the moments, as well as the as deter- 
mination, follow closely the analysis by the OPAL experi- 
ment in the complete LEP energy range of 91-209 GeV [7]. 
This work supplements our previous analyses on jet pro- 
duction rates, determinations of as and four jet produc- 
tion, using JADE and OPAL data [3,8,9]. 

The outline of the paper is as follows. In Sect. [21 we 
present the observables used in the analysis and describe 
the perturbative QCD predictions. In Sect. [3] the analy- 



sis procedure is explained in detail. Sect. [4] contains the 
discussion of the systematic checks which are performed 
and the resulting systematic errors. We collect the results 
and describe the determination of as in Sect. [5l and we 
summarize in Sect. [B] 



2 Observables 

Event shape variables are a convenient way to characterise 
properties of hadronic events by the distribution of parti- 
cle momenta. For the definition of the variables we refer 
to [7]. The event shapes considered here are Thrust T, 
C-parameter C, Heavy Jet Mass Mh, jet broadening vari- 
ables Bt and Byf, and the transition value between 2 and 

3 jets in the Durham jet scheme, 2/^3. The as determi- 
nation in [7] is based on distributions and moments of 
these variables. Their theoretical description is currently 
the most advanced [10-12]. Further, we measure moments 
of Thrust major Tmaj, Thrust minor Tmin, Oblateness O, 
Sphericity S, Light Jet Mass Ml, and Narrow Jet Broad- 
ening Bjij. Moments of these variables and variances of 
all measured event shapes will be made available in the 
HEPDATA databaseQ. 

Generic event shape variables y are constructed such 
that spherical and multi-jet events yield large values of y, 
while two narrow back-to-back jets yield ?/ ~ 0. Thrust T 



^ |http : //durpdg . dur . ac . uk/HEPDATA/ | 
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is an exception to this rule. By using y = 1 — T instead 
the condition is fulfilled for all event shapes. 

The Tith, n = 1,2,..., moment of the distribution of 
the event shape variable y is defined by 

j-y^a^ 1 dcr 

Jo cr dy 

where ymax is the kinematically allowed upper limit of the 
variable y. 

Predictions have been calculated for the moments of 
event shapes. Their evolution with cm. energy allows di- 
rect tests of the predicted energy evolution of the strong 
coupling as . Furthermore we determine as (-^z" ) by evolv- 
ing our measurements to the energy scale given by the 
mass of the Z° boson. The theoretical calculations involve 
a integration over full phase space, which implies that 
comparison with data always probes all of the available 
phase space. This is in contrast to QCD predictions for dis- 
tributions; these are commonly only compared with data- 
e.g. in order to measure as-in restricted regions, where 
the theory is well defined and describes the data well, see 
e.g. [2]. Comparisons of QCD predictions for moments of 
event shape distributions with data are thus complemen- 
tary to tests of the theory using distributions. 

Uncertainties in the NNLO predictions for event shape 
distributions in the two-jet region [11,12] prevent the re- 
liable calculation of moments to NNLO at present, and 
therefore we compare with NLO predictions only. The 
QCD prediction of (y") at parton level, in next-to-leading 
order (NLO) perturbation theory, and with as = as/(27r), 
is 

^^„^part.theo ^ + (g^^ „ 2^„) . (2) 

The values of the coefficienttQ An and Bn can be obtained 
by numerical integration of the QCD matrix elements us- 
ing the program EVENT2 [13]. 

The coupling as and the ag coefficient depend on the 
renormalisation scale /i [14]. For the sake of clarity the 
renormalisation scale factor is defined as = [xjsfs^ so 
setting to one implies that the renormalisation scale is 
y/s. A truncated fixed order QCD calculation such as ([2]) 
will then depend on . The renormalisation scale depen- 
dence is implemented by the replacement S„ — > Bn + 
/3o ln(a;^)y^„ where /3o = 11 — ffif is the leading order /3- 
function coefficient of the renormalisation group equation 
and nt = 5 is the number of active quark ftavours. 



3 Analysis procedure 

3.1 The JADE detector 

The JADE detector is described in detail in ref. [1]. En- 
ergy measurement by the electromagnetic calorimeter and 
the reconstruction of charged particle tracks in the central 

^ The a| coefficient is written as Bn — 2An because the QCD 
calculations are normalized to the Born cross section ctq, while 
the data are normalized to the total hadronic cross section, 
cTtot = cro(l + 2as) in LO. 



track detector are the main ingredients for this analysis. 
The central jet chamber was positioned in an axial mag- 
netic field of 0.48 T provided by a solenoidal magnet0 The 
magnet coil was surrounded by the lead glass calorimeter, 
which measured electromagnetic energy and consisted of 
a barrel and two endcap sections. 

3.2 Data samples 

In this analysis we are using data samples identical to the 
samples used in [1-4,8,15], collected by the JADE exper- 
iment between 1979 and 1986; they correspond to a total 
integrated luminosity of ca. 195 pb~^. Table[T]contains the 
breakdown of the data samples-data taking period, energy 
range, mean centre-of-mass energy, integrated luminosity 
and the number of selected hadronic events. 



Table 1. Year of data taking, energy range, integrated lu- 
minosity, average centre-of-mass energy and the numbers of 
selected data events for each data sample 



year 




range of 
^ in GeV 


^/s mean 
in GeV 


luminosity 
(pb-i) 


selected 
events 


1981 




13.0-15.0 


14.0 


1.46 


1783 


1981 




21.0-23.0 


22.0 


2.41 


1403 


1981- 


-1982 


33.8-36.0 


34.6 


61.7 


14313 


1986 




34.0-36.0 


35.0 


92.3 


20876 


1985 




37.3-39.3 


38.3 


8.28 


1585 


1984- 


-1985 


43.4-46.4 


43.8 


28.8 


4376 



3.3 Monte Carlo samples 

To correct the data for experimental effects and back- 
grounds we use samples of MC simulated events. Using 
PYTHIA 5.7 [16] we simulate the process e"'"e~ — > hadrons. 
For systematic checks we use corresponding samples ob- 
tained by simulation with HERWIG 5.9 [17]. We process 
the MC samples generated at each energy point through a 
full simulation of the JADE detector [18-20], summarized 
in [15]; and we reconstruct them in essentially the same 
way as the data. 

Using the parton shower models PYTHIA 6.158, HER- 
WIG 6.2 [21] and ARIADNE 4.11 [22] we employ in ad- 
dition large samples of MC events without detector sim- 
ulation, in order to compare with the corrected data. For 
the purpose of comparison with the data, the MC events 
include the effects of hadronisation, i.e. the transition of 
partons into hadrons. All used major versions of the mod- 
els were adjusted to LEP 1 data by the OPAL collabora- 
tion [23,24], so we expect comparable results from them. 

In the JADE right-handed coordinate system the +x axis 
pointed towards the centre of the PETRA ring, the y axis 
pointed upwards and the z axis pointed in the direction of the 
positron beam. The polar angle and the azimuthal angle (f) 
were defined with respect to z and x, respectively, while r was 
the distance from the z-axis. 
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3.4 Selection of events 

The selection-identical to the one used in [8]-aims at se- 
lecting hadronic events in the JADE data excluding events 
with much energy lost by initial state radiation (ISR). The 
rejected background consists to a large degree of two pho- 
ton events. It uses cuts on event multiplicity, on visible 
energy and longitudinal momentum balance. The cuts are 
documented in [2,25-27]. 

So called good tracks and calorimeter clusters are iden- 
tified by appropriate criteria [8]. Double counting of en- 
ergy from charged tracks and calorimeter clusters is avoi- 
ded by subtracting the estimated contribution of a charged 
track from the associated cluster energy. 

The number of selected events for each energy point is 
given in table [TJ 



3.5 Corrections to the Data 




10 20 30 40 10 20 30 40 

Vs"{GeV) Vs(GeV) 



The data are corrected further for the effects of limited 
detector acceptance and resolution, and residual ISR fol- 
lowing [8]. All selected charged tracks, as well as the elec- 
tromagnetic calorimeter clusters remaining after the cor- 
rection for double counting of energy as described above, 
are used in the evaluation of the event shape moments. 
The values of the moments after the application of all se- 
lection cuts are said to be at the detector level. 

As the QCD predictions are calculated for massless 
quarks we have to correct our data for the presence of 
events originating from bb final states. Especially at low 
y/s the large mass of the b quarks and of the subsequently 
produced and decaying B hadrons will influence the values 
of the event shape variables. Therefore in the JADE anal- 
ysis events from the process e^e~ bb (approximately 
1/11 of the hadronic events) are considered as background. 

For the determination of the moments we calculate the 
sums Vi'da.ts. i^*^^ moment order n — 1, . . . , 5) where i 
counts all selected events. The expected contribution of bb 
background events X^i^rbb' ^ estimated by PYTHIA, is 
subtracted from the observed sum 2/"data- ^ multi- 
plicative correction we then account for the effects of de- 
tector imperfections and of residual ISR and two photon 
background. 

Two sets of sums ^ ■ are calculated from MC sim- 
ulated signal events. At detector level, MC events are 
treated identically to the data. The hadron level set is 
computed using the true momenta of the stable parti- 
cles in the eventQ, and uses only events where \fs' , the 
cm. energy of the event, reduced due to ISR, satisfies 
Vs - < 0.15 GeV. The ratio of the MC hadron level 
moment over the MC detector level moment is applied as a 
detector correction factor for the data; the corrected sums 
are normalized by the corrected total number of selected 



* All charged and neutral particles with a lifetime larger 
than 3 x 10"^" s are considered stable. 



Fig. 1. Detector correction factors Cdct as calculated using 
the PYTHIA MC model (see text for details). Line types cor- 
respond to moment order as shown in top left figure 

events A'tot yielding the final values of (j/"), 

, „vhad.MC / \ 

(3) 

The corrected total number of events is calculated from 
the number of selected events in the data in the same way 
as for the moments. 

There is some disagreement between the detector cor- 
rections calculated using PYTHIA or HERWIG at low 

while at larger ^/s the correction factors agree well for 
most observables. The difference in detector corrections 
will be evaluated as an experimental systematic uncer- 
tainty, see Sect. |4l The detector correction factors 
CdGt= as determined using PY- 

THIA are shown in Fig. [TJ 



4 Systematic uncertainties 

Several contributions to the experimental uncertainties 
are estimated by repeating the analysis with varied track 
or event selection cuts or varied procedures as in [8]. For 
each systematic variation the value of the event shape mo- 
ment or of as is determined and then compared to the 
default value. The experimental systematic uncertainty 
quoted is the result of adding in quadrature all contri- 
butions. In the fits of the QCD predictions to the data 
two further systematic uncertainties are evaluated: 

- Using HERWIG 6.2 and ARIADNE 4.11 instead of 
PYTHIA 6.158 we assess the uncertainties associated 
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with the hadronisation correction (Sect. [5T2|l . The hadro- 
nisation systematic uncertainty is defined by the larger 
change in as resulting from these alternatives. 
By varying the renormalisation scale factor a;^ we as- 
sess the theoretical uncertainty associated with miss- 
ing higher order terms in the theoretical prediction. 
The renormahsation scale factor is set to 2.0 and 
0.5. The theoretical systematic uncertainty is defined 
by the larger deviation from the default value. 



5 Results 

5.1 Values of event shape moments 

The first five moments of the six event shape variables 
after subtraction of bb background and correction for de- 
tector effects measured by JADE are listed in Tables [5] 
and [3] and shown in Figs. [2] and [3l Superimposed we show 
the moments predicted by the PYTHIA, HERWIG and 
ARIADNE MC models tuned by OPAL to LEP 1 data. 
The moments become smaller by approximately one or- 
der of magnitude with each increasing moment order; the 
higher moments are more strongly suppressed with centre- 
of-mass energy. Statistical and experimental systematic 
uncertainties strongly increase with moment order. 

In order to make a clearer comparison between data 
and models the lower plots in Figs. [2] and [3] show the dif- 
ferences between data and each model divided by the com- 
bined statistical and experimental error for ^/s — 14 and 
35 GeV. The three models are seen to describe the data 
fairly well; PYTHIA and ARIADNE are found to agree 
better with the data than HERWIG. 



5.2 Determination of cts 

In order to measure the strong coupling as, we fit the 
QCD predictions to the corrected moment values (y"), 
i.e. to the data shown in Tables O and [31 The theoreti- 
cal predictions using the 0(ag) calculation described in 
Sect. [2] provide values at the parton level. It is necessary 
to correct for hadronisation effects in order to compare the 
theory with the hadron level data. Therefore the moments 
are calculated at hadron and parton level using large sam- 
ples of PYTHIA 6.158 events and, as a cross check, sam- 
ples obtained by simulation with HERWIG 6.2 and ARI- 
ADNE 4.11. Parton level is the stage at the end of the 
parton shower in the simulation of an hadronic event. In 
order to correct for hadronisation the data points are then 



multiplied by the ratio C^""^ = (y")Part,MC^^y„^ 
the parton over hadron level moments; ^y"^P'*'^' 



had,MC 



of 



C 



had 



had 



The models use cuts on quantities like e.g. the invari- 
ant mass between partons in order to regulate divergencies 
in the predictions for the parton shower evolution. As a 
consequence in some events no parton shower is simulated 
and the original quark-antiquark pair enters the hadroni- 
sation stage of the model directly. This leads to a bias in 



(d-T)") 



n=3 
n=4 
n=5 



40 



20 30 40 

*(GeV) 




30 40 

Vs(GeV) 



Fig. 4. Hadronisation correction factors C^'^'^ as calculated 
using the MC model PYTHIA 6.158 (see text for details). Line 
types correspond to moment order as shown in top left figure 



the calculation of moments at the parton level, since y = 
in this case for all observables considered here (?/^ cannot 
be calculated in this case). In order to avoid this bias we 
exclude in the simulation at the parton level events with- 
out gluon radiation, as in [28]. After this exclusion, the 
^/s evolution of the moments follows the QCD prediction; 
the change of the prediction is comparable in size with 
the differences between employed MG generators. At the 
hadron and detector level all events are used. 

The hadronisation correction factors C^^'^ as obtained 
from PYTHIA 6.158 are shown in Fig.Hl We find that the 
hadronisation correction factors can be as large as two 
at low y/s. For larger y/s the hadronisation corrections 
decrease as expected. 

The models PYTHIA 6.158, HERWIG 6.2 and ARI- 
ADNE 4.11 do not agree well for moments based on _Bw, 

and Afn at low ^/s. The differences between the models 
are studied as a systematic uncertainty in the fits. 

A value for each moment (?/") is calculated using 
the formula 



X 



, part 



1 part,thoo\2 / 2 
/i ) l^i ■ 



(4) 



where i counts the energy points, <Ti denotes the statistical 
errors and ^y")P*^''*'*'^''° ig taken from 

The value is minimized with respect to as(Mzo) 
for each moment n separately. The statistical uncertainty 
is found by varying the minimum value to Xmin + 

1. The evolution of as from Mip to cm. energy {^/s)i 
is implemented in the fit in two-loop precision [29]. The 
renormalisation scale factor x^, as discussed in Sect. [U is 
set to 1. 
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Fig. 2. First five moments of 1 — T, C and Bt at hadron level compared with predictions based on PYTHIA 6.158, HERWIG 6.2 
and ARIADNE 4.11 MC events. The inner error bars-where visible-show the statistical errors, the outer bars show the total 
errors. Where no error bar is visible, the total error is smaller than the point size. The lower panels show the differences between 
data and MC at = 14 and 35 GeV, divided by the total error 
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Fig. 3. First five moments of Bw, J/23 and Mh at hadron level compared with predictions based on PYTHIA 6.158, HERWIG 6.2 
and ARIADNE 4.11 MC events. The inner error bars-where visible-show the statistical errors, the outer bars show the total 
errors. Where no error bar is visible, the total error is smaller than the point size. The lower panels show the differences between 
data and MC at = 14 and 35 GeV, divided by the total error 
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5.3 Fits of JADE data 

Data and NLO prediction are compared, and this is re- 
peated for every systematic variation. The results are 
shown in Fig. [5] and listed in Table S) Figure [5] also con- 
tains the combination of the fit results discussed below. 
The values of x^/d.o.f . are in the order of 1-10, the fitted 
predictions-including the energy evolution of as-are con- 
sistent with the data. The fit to (M^) does not converge 
and therefore no result is shown □ 
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Fig. 5. Measurements of aslPI^o) using fits to moments of 
six event shape variables at PETRA energies. The inner error 
bars-where visible-show the statistical errors, the outer bars 
show the total errors. The dotted line indicates the weighted 
average described in Subsect. 15.51 the shaded band shows its 
error. Only the measurements indicated by solid symbols are 
used for this purpose 

The fitted values of as (Mzo ) increase steeply with the 
order n of the moment used, for ((1 — T)"), (C") and 
(_B^). This effect is less pronounced and systematic for 
(5w>> ((ySi)") and (Mg). In Fig.[9]we show the ratio K = 
Bn/An of NLO and LO coefficients for the six observables 
used in our fits to investigate the origin of this behaviour. 
Steeply increasing values of as{Mzo) with moment order 
n for ((1 — r)"), (C") and (BJ^) and increasing values of 
K with n are clearly correlated. There is also a correlation 
with the rather large scale uncertainties in the respective 

^ Equation [2] precludes a real solution as , if Bn — '2 An < 
-Al/4:{y"). For (Mh) this relation is fulfilled in the whole 
energy range 14-207 GeV, see Tables [2] and |3] and [7]. The 
NLO coefficient is negative in the case of (Bw)> too. This ob- 
servable gives the maximum value of x^/d.o.f .=98.5/5, further 
problems in the determination of as using (Bw) show up in 
Subsect. [531 



fits. The other observables {B^), ((ygj)"> and (Mg) have 
more stable results for ag {Mio ) and correspondingly fairly 
constant values of K. The reason that the fit of (Af^) does 
not converge is the large and negative value of K. 

5.4 Combined fits of JADE and OPAL data 

For the most significant results we supplement the JADE 
data with the analogous OPAL data [7], covering the en- 
ergy range of 91 to 209 GeV. 

The JADE and OPAL detectors are very similar, both 
in construction and in the values of many detector param- 
eters. The combined use of the JADE and OPAL data can 
therefore be expected to lead to consistent measurements, 
with small systematic differences. Our analysis procedure 
is therefore constructed to be similar to the one used in 
the OPAL analysis [7] , in particular in the estimate of the 
systematic errors. 

The central values and statistical errors of the com- 
bined fits are found employing the calculation ([4]) as 
aboveQ However, the systematic uncertainties in this case 
cannot be found by simple repetitions of the fits, as the 
JADE and OPAL systematic variations are not identical. 

The systematic uncertainties are correlated between 
different energy points, and including general correlations, 
the calculation shown in ([4]) has to be generalised to 
[30] 

x'-E ((2/">r - (y">r"''=°) • 

where the Vij are the covariances of the n-th moment 
at the energy points i and j. They have the form Vij = 
Sij + Eij , with statistical covariances Sij and experimen- 
tal systematic covariances Eij. The matrix Sij is diagonal, 
Sii = CTg^at i 1 as data of different energy points are inde- 
pendent. The experimental systematic covariances Eij are 
only partly known: 

— The diagonal entries are given by 

denoting by crcxp..i the experimental uncertainty at en- 
ergy point i. 

— The non diagonal entries can only follow from plau- 
sible assurnptions. We employ the minimum overlap 
assumptioT^ 

E,j = Min{CT^^p , cr^xp.,j} > (6) 

® For this reason systematic differences between the two 
experiments contribute to the sometimes high values; in 
Figs. Eland [7] the experimental uncertainties are indicated sep- 
arately. 

^ Fitting the low energy JADE points (14, 22 GeV) this 
assumption results [31] in a more accurate and more con- 
servative error estimation than the full overlap assumption 
E,j = Max{(T^xp.,i , employed in [3]. 
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Fig. 6. Fits of the NLO predictions to JADE and OPAL mea- 
surements of moments of C at parton leveL The solid lines 
show the -^s evolution of the NLO prediction based on the fit- 
ted value of asiM^o). The inner error bars-where visible-show 
the statistical errors used in the fit, the outer bars show the 
total errors. Where no error bar is visible, the total error is 
smaller than the point size 



Fig. 7. Fits of the NLO predictions to JADE and OPAL mea- 
surements of moments of Bw at parton level. The solid lines 
show the -^s evolution of the NLO prediction based on the fit- 
ted value of as {M^o ) . The inner error bars-where visible-show 
the statistical errors used in the fit, the outer bars show the 
total errors. Where no error bar is visible, the total error is 
smaller than the point size. Problems of the NLO prediction 
at low -v/s are discussed in the text 



The total errors are found by fits employing the 
calculation ([5]). We use the relative experimental uncer- 
tainties to determine the experimental uncertainties of the 
central values from the fits without correlations. 

Figures [S] and [7] show the comparison of data points 
and predictions for the moments of the C-parameter and 
the wide jet broadening _Bw The predictions for (C") 
are seen to be in good agreement with the data and sig- 
nificantly confirm the QCD prediction of the energy de- 
pendence of as(-\/s), also known as asymptotic freedom. 
The prediction slightly overshoots the higher moments of 
1 — T, C and Bt at 14GeV, and undershoots the mo- 
ments of i?W: -^^H, and sometimes 2/^3. At low y/s the 
NLO predictions of the i?w, 2/23 and Afn distributions are 
(unphysically) negative in a large range of the two jet re- 
gion [15]. Therefore the NLO prediction for the moments 
is likely to be incomplete and too low to provide a sat- 
isfactory description of the data at low cm. energies. In 
the case of {B^) the coefficient is even negative, and 
we do not show or use this fit. The results are listed in 
Tableland shown in Fig.[8l 



5.5 Combination of o:s determinations 

To make full use of the data, we combine the measure- 
ments of as{Mzo) from the various moments and event 
shapes and determine a single value. An extensive study 



was done by the LEP QCD working group on this prob- 
lem [6,7,32-34], and their procedure is adopted here. 

A weighted mean of the ag (Mzo ) measurements is cal- 
culated which minimizes the formed from the measure- 
ments and the combined value. This mean value, as (-^z" ) i 
is given by 

as(Afzo) = ^Wias,i with Wi = , 

(7) 

where the measured values of asC-^z") are denoted as,i, 
their covariance matrix V\ and the individual results are 
counted by «, j and k. The presence of highly correlated 
and dominant systematic errors makes a reliable estimate 
of V difficult. Undesirable features (such as negative 
weights) can be caused by small uncertainties in the esti- 
mation of these correlations. Therefore only experimental 
systematic errors-assumed to be partially correlated by 
minimum overlap as Vl^ = min(cr^xp j, fj^^p j)~and statisti- 
cal correlations are taken to contribute to the off-diagonal 
elements of the covariance matrix. The statistical correla- 
tions are determined using MC simulations at the parton 
levelH The diagonal elements are calculated from all error 
contributions-statistical, experimental, hadronisation and 

* The result is identical if the correlations are determined us- 
ing PYTHIA, HERWIG or ARIADNE at 14.0.. .43.8 GeV, or 
determined at hadron level instead of parton level. The corre- 
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Fig. 8. Measurements of qs(Mzo) using fits to moments of six 
event shape variables at PETRA and LEP energies. The inner 
error bars-where visible-show the statistical errors, the outer 
bars show the total errors. The experimental systematic un- 
certainties are estimated by the minimum overlap assumption. 
The dotted line indicates the weighted average described in the 
text, the shaded band shows its error. Only the measurements 
indicated by solid symbols are used for this purpose 

theory uncertainties. Using the weights derived from the 
covariance matrix V' the theory uncertainties are com- 
puted by analogously combining the as{Mzo) values from 
setting = 2.0 or = 0.5, and the hadronisation un- 
certainties by combining the results obtained with the al- 
ternative hadronisation models. 

To select observables with an apparently converging 
perturbative prediction, we consider [7] only those results 
for which the NLO term in equation ([2]) is less than half 
the corresponding LO term (i.e. iKas/^Trl < 0.5 or \K\ < 
25), namely (1 - T), (C), (St), {B^) and {{ygr), n = 
1, . . . , 5; and (M^), n = 2, . . . , 5. These are results from 
17 observables in total; or 16 observables from JADE and 
OPAL, excluding (B^). The K values are shown in Fig.O 

Using only JADE data, the result of the combination 

is 

as(Mzo) = 0.1287 ± 0.0007(stat.) ± 0.0011(exp.) 
±0.0022(had.) ± 0.0075(theo.) , 

and is shown in Fig. [5] Combining JADE and OPAL mea- 
surements, the result is 

as(Mzo) = 0.1262 ±0.0006(stat.)±0.0010(exp.) 
±0.0007(had.) ± 0.0064(theo.) , 

lation values are cited in [31]; at 14 GeV and parton level they 
vary between 29% and 99% and are larger than 50% mostly. 
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Fig. 9. The ratio K = Bn/A„ of NLO and LO coeflicients for 
the first five moments of the six event shape variables used in 
the determination of as, see also [7] 

and is shown in Fig. [S] Both values are above, but still 
consistent with the world average of as(Mzo) = 0.1189 ± 
0.0010 [35]. It has been observed previously in compar- 
isons of event shape distributions with NLO [36] or 
NNLO [37] QCD predictions with a;^ = 1 that fitted val- 
ues of as{Mzo) tend to be large compared to the world 
average. 

To enable comparison with earlier and more specific 
analyses [38] we combine the JADE fit results from only 
the firsiEl moments (1-T), (C), (Bt), {Bw), (yFa) and 
(A/y). This yields a value of 

as(Afzo) = 0.1243 ± 0.0001(stat.) ± 0.0009(exp.) 
±0.0010(had.) ± 0.0070(theo.) . 

The slightly smaller error in this determination of as re- 
flects the fact that the lower order moments are less sensi- 
tive to the multijet region of the event shape distributions. 
This leads to a smaller statistical and theoretical uncer- 
tainty. In all three measurements the scale uncertainty is 
dominant. 



6 Summary 

In this paper we present measurements of moments of 
event shape distributions at centre-of-mass energies be- 
tween 14 and 44 GeV using data of the JADE exper- 
iment. The predictions of the PYTHIA, HERWIG and 
ARIADNE MC models tuned by OPAL to LEP 1 data 

^ Because of the problems with the NLO description of 
(M^y'^^, {M|) is often regarded as the first moment of Mh. 
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are found to be in reasonable agreement with the mea- 
sured moments. 

From fits of 0{a'g) predictions to selected event shape 
moments corrected for experimental and hadronisation ef- 
fects we have determined the strong coupling to be 
as(Mzo) = 0.1287±0.0079 (total error) using only JADE 
data, and as(Mzo) = 0.1262 ± 0.0065 (total error) us- 
ing combined JADE and OPAL data. Fits to moments 
of Mh, Bw and j/23 return large values of x^/d-o.f.; the 
higher moments, in particular of the 1— T, C and Bt event 
shape variables, lead to systematically enlarged values of 
as- Results where ag is steeply rising with moment order 
are strongly correlated with the relative size of the Oq co- 
efficient and thus are most likely affected by deficiencies 
of the NLO prediction. 

The JADE experiment assesses an interesting energy 
range for the perturbative analysis since the energy evo- 
lution of the strong coupling is more pronounced at low 
energies. 
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Table 2. Moments of the 1 - T, C, Bt, Bw, 2/23 and Mh distributions measured by JADE at 14.0, 22.0 and 34.6 GeV. The 
first uncertainty is statistical, while the second is systematic 



((1 - T)") at 14.0 GeV 



((1 - T)") at 22.0 GeV 



((1 - T)") at 34.6 GeV 



1 


(1.405 ± 0.022 


± 0.050) 


10" 


-1 


(1.123 ± 0.021 ± 0.028) 


10" 


-1 


2 


(2.38 ± 0.08 


± 0.17 ) 


10" 


-2 


(1.700 ± 0.068 ± 0.086) 


10" 


-2 


3 


(4.68 ± 0.28 


± 0.54 ) 


10" 


-3 


(3.31 ± 0.21 ± 0.26 ) 


10" 


-3 


4 


(1.04 ± 0.09 


± 0.17 ) 


10" 


-3 


(7.79 ± 0.69 ± 0.84 ) 


10" 


-4 


5 


(2.55 ± 0.33 


± 0.56 ) 


10" 


-4 


(2.08 ± 0.24 ± 0.29 ) 


10" 


-4 



(8.99 ± 0.07 ± 0.13 )-10-^ 
(1.192 ± 0.020 ± 0.024)-10"^ 
(2.151 ± 0.057 ± 0.052)-10-^ 
(4.77 ± 0.17 ± 0.16 )-10-'' 
(1.202 ± 0.056 ± 0.061)-10"'' 



(C") at 14.0 GeV 



(C") at 22.0 GcV 



(C") at 34.6 GoV 



1 (5.22 ± 0.05 ± 0.13 )-10-^ 

2 (3.00 ± 0.06 ± 0.14 )-10-^ 

3 (1.85 ± 0.06 ± 0.13 )-10-^ 

4 (1.22 ± 0.05 ± 0.11 )-10-^ 

5 (8.39 ± 0.45 ± 0.96 j-lO'^ 



(4.280 ± 0.057 ± 0.077)- 10"^ 
(2.152 ± 0.056 ± 0.075)-10-^ 
(1.235 ± 0.047 ± 0.064)-10-^ 
(7.86 ± 0.39 ± 0.54 j-lO"^ 
(5.41 ± 0.33 ± 0.46 j-lQ-^ 



(3.512 ± 0.020 ± 0.043)-10-i 
(1.561 ± 0.018 ± 0.030)-10-^ 
(8.31 ± 0.14 ± 0.19 )-10-^ 
(5.04 ± 0.11 ± 0.13 )-10-^ 
(3.335 ± 0.088 ± 0.099)-10-2 



(S?) at 14.0 GeV 



(B?) at 22.0 GeV 



(B?) at 34.6 GeV 



n 



(1.918 ± 0.017 ± 0.038)40"^ 
•10"^ 
•10"^ 

■10"=* 
■10"'' 



(3.94 ± 0.07 ± 0.16 

(8.61 ± 0.26 ± 0.54 

(1.99 ± 0.09 ± 0.17 

(4.84 ± 0.28 ± 0.55 



(1.627 ± 0.018 ± 0.021)-10"^ 
(2.963 ± 0.067 ± 0.077)-10'^ 
(6.01 ± 0.21 ± 0.24 )-10-^ 
(1.344 ± 0.065 ± 0.072)-10-^ 
(3.26 ± 0.20 ± 0.22 )-10-* 



(1.372 ± 0.006 ± 0.011)-10"^ 
(2.202 ± 0.021 ± 0.033)-10~^ 
(4.082 ± 0.062 ± 0.083)-10-^ 
(8.56 ± 0.18 ± 0.22 j-lO""* 
(1.978 ± 0.053 ± 0.062)-10-* 



(B^) at 14.0 GeV 



(B^) at 22.0 GeV 



(B^) at 34.6 GeV 



1 (1.166 ± 0.011 ± 0.019)-10"^ 

2 (1.482 ± 0.031 ± 0.048)-10"^ 

3 (2.045 ± 0.070 ± 0.098)-10-^ 

4 (3.04 ± 0.15 ± 0.19 )-10-'* 

5 (4.81 ± 0.33 ± 0.38 )-10-^ 



n 



(1.000 ± 0.012 ± 0.014)-10"^ 
(1.151 ± 0.030 ± 0.033)-10"^ 
(1.525 ± 0.065 ± 0.068)-10"^ 
(2.30 ± 0.14 ± 0.15 )-10-* 
(3.84 ± 0.32 ± 0.36 j-lO"^ 



(8.720 ± 0.047 ± 0.087)-10-^ 
(9.42 ± 0.11 ± 0.20 )-10"^ 
(1.238 ± 0.023 ± 0.042)-10-3 
(1.897 ± 0.050 ± 0.085)-10-'' 
(3.25 ± 0.11 ± 0.18 )-10-^ 



{{y&r) at 14.0 GeV 



((2/?3)") at 22.0 GeV 



((j/?3)") at 34.6 GeV 



1 


(3.54 ± 0.12 ± 0.19 ) 


10" 


-2 


(2.89 


± 0.12 


± 0.10 ) 


10" 


-2 


(2.408 ± 0.042 ± 0.041)-10" 


-2 


2 


(2.55 ± 0.22 ± 0.29 ) 


10" 


-3 


(2.55 


± 0.23 


± 0.18 ) 


10" 


-3 


(2.173 ± 0.081 ± 0.049)-10" 


-3 


3 


(3.01 ± 0.48 ± 0.48 ) 


10" 


-4 


(3.93 


± 0.55 


± 0.37 ) 


10" 


-4 


(3.35 ± 0.19 ± 0.13 )-10" 


-4 


4 


(4.7 ± 1.2 ± 1.0 ) 


10" 


-5 


(7.6 


± 1.4 


± 0.8 ) 


10" 


-5 


(6.42 ± 0.48 ± 0.38 )-10" 


-5 


5 


(8.8 ± 3.0 ± 2.5 ) 


10" 


-6 


(1.66 


± 0.36 


± 0.21 ) 


10" 


-5 


(1.38 ± 0.13 ± 0.11 )-10" 


-6 


n 


(MS) at 14.0 GeV 








(MS) at 22.0 GeV 






(MS) at 34.6 GeV 




1 


(3.207 ± 0.024 ± 0.049) 


10" 


-1 


(2.832 ± 0.026 


± 0.036) 


10" 


-1 


(2.522 ± 0.010 ± 0.024)-10" 


-1 


2 


(1.074 ± 0.017 ± 0.033) 


10" 


-1 


(8.55 


± 0.16 


± 0.21 ) 


10" 


-2 


(6.979 ± 0.057 ± 0.095)-10" 


-2 


3 


(3.74 ± 0.09 ± 0.17 ) 


10" 


-2 


(2.74 


± 0.08 


± 0.10 ) 


10" 


-2 


(2.114 ± 0.028 ± 0.030)-10" 


-2 


4 


(1.348 ± 0.048 ± 0.085) 


10" 


-2 


(9.36 


± 0.40 


± 0.49 ) 


10" 


-3 


(6.98 ± 0.13 ± 0.10 )-10" 


-3 


5 


(5.02 ± 0.24 ± 0.41 ) 


10" 


-3 


(3.38 


± 0.20 


± 0.23 ) 


10" 


-3 


(2.485 ± 0.062 ± 0.051)-10" 


-3 
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Table 3. Moments of the 1 - T, C, Bt, Bw, 2/?3 and Mh distributions measured by JADE at 35.0, 38.3 and 43.8 GeV. The 
first uncertainty is statistical, while the second is systematic 



((1 - T)") at 35.0 GeV 



((1 - T)") at 38.3 GeV 



((1 - T)") at 43.8 GeV 



1 (9.22 ± 0.07 ± 0.18 

2 (1.260 ± 0.019 ± 0.045) 

3 (2.34 ± 0.06 ± 0.11 

4 (5.36 ± 0.17 ± 0.29 

5 (1.394 ± 0.058 ± 0.084) 




(8.07 ± 0.12 ± 0.10 )-10-^ 
(1.032 ± 0.032 ± 0.024)-10-^ 
(1.867 ± 0.093 ± 0.069)-10-^ 
(4.22 ± 0.29 ± 0.23 j-lO"" 
(1.099 ± 0.098 ± 0.084)-10-'' 



(C") at 35.0 GcV 



(C") at 38.3 GcV 



(C") at 43.8 GoV 



1 (3.582 ± 0.019 ± 0.057)-10"^ 

2 (1.620 ± 0.017 ± 0.047)-10-^ 

3 (8.76 ± 0.13 ± 0.34 )-10-^ 

4 (5.38 ± 0.10 ± 0.25 j-lO"^ 

5 (3.60 ± 0.08 ± 0.19 j-lQ-^ 



(3.486 ± 0.056 ± 0.065)-10-^ 
(1.587 ± 0.048 ± 0.052)-10"^ 
(8.76 ± 0.38 ± 0.39 )-10-^ 
(5.49 ± 0.30 ± 0.30 j-lO"^ 
(3.74 ± 0.25 ± 0.24 j-lO'^ 



(3.178 ± 0.034 ± 0.032)- 10^1 
(1.347 ± 0.028 ± 0.025)-10-^ 
(7.05 ± 0.22 ± 0.18 )-10-^ 
(4.25 ± 0.17 ± 0.14 )-10-^ 
(2.81 ± 0.14 ± 0.11 )-10-2 



(S?) at 35.0 GeV 



(B?) at 38.3 GeV 



(B?) at 43.8 GeV 



1 (1.395 ± 0.006 ± 0.016) 

2 (2.277 ± 0.020 ± 0.054). 

3 (4.30 ± 0.06 ± 0.15 

4 (9.20 ± 0.17 ± 0.42 

5 (2.17 ± 0.05 ± 0.12 



n 



10" 


-1 


10" 


-2 


10" 


-3 


10" 


-4 


10" 


-4 



(1.364 ± 0.018 ± 0.019) 
(2.229 ± 0.059 ± 0.064) 
(4.27 ± 0.17 ± 0.18 ) 
(9.30 ± 0.51 ± 0.54 ) 
(2.23 ± 0.15 ± 0.16 ) 



10" 


-1 


10" 


-2 


10" 


-3 


10" 


-4 


10" 


-4 



(1.260 ± 0.011 ± 0.010)-10"^ 
(1.920 ± 0.035 ± 0.031)-10-2 
(3.480 ± 0.098 ± 0.087)-10-^ 
(7.26 ± 0.28 ± 0.26 j-lO""* 
(1.681 ± 0.084 ± 0.080)-10-* 



(B^) at 35.0 GeV 



(B^) at 38.3 GeV 



(B^) at 43.8 GeV 



n 



1 (8.90 ± 0.04 ± 0.13 )-10-'' 

2 (9.83 ± 0.10 ± 0.31 ^lO"'* 

3 (1.323 ± 0.022 ± 0.061)-10-=^ 

4 (2.08 ± 0.05 ± 0.12 )-10-'* 

5 (3.65 ± 0.11 ± 0.24 )-10-^ 



(8.76 ± 0.13 ± 0.17 )-10-^ 
(9.72 ± 0.29 ± 0.40 )-10"''* 
(1.319 ± 0.062 ± 0.087)-10-^ 
(2.07 ± 0.13 ± 0.20 )-10-* 
(3.60 ± 0.30 ± 0.50 )-10-^ 



(8.185 ± 0.082 ± 0.097)-10-^ 
(8.73 ± 0.19 ± 0.22 )-10"^ 
(1.171 ± 0.040 ± 0.046)-10-3 
(1.859 ± 0.089 ± O.lOoj-lO-" 
(3.31 ± 0.21 ± 0.23 )-10-^ 



((2/23)") at 35.0 GeV 



((2/?3)") at 38.3 GeV 



((j/?3)") at 43.8 GeV 



1 


(2.551 ± 0.040 ± 0.058) 


10" 


-2 


(2.66 


± 0.12 


± 0.16 ) 


10" 


-2 


(2.269 ± 0.071 


=t 0.068)-10-^ 


2 


(2.395 ± 0.080 ± 0.071) 


10" 


-3 


(2.62 


± 0.23 


± 0.28 ) 


10" 


-3 


(2.10 ± 0.14 


± 0.11 )-10"^ 


3 


(3.77 ± 0.19 ± 0.23 ) 


10" 


-4 


(4.04 


± 0.52 


± 0.61 ) 


10" 


-4 


(3.24 ± 0.32 


± 0.28 )-10"'' 


4 


(7.31 ± 0.48 ± 0.71 ) 


10" 


-5 


(7.5 


± 1.3 


± 1.6 ) 


10" 


-5 


(6.19 ± 0.82 


± 0.77 )-10-^ 


5 


(1.58 ± 0.13 ± 0.20 ) 


10" 


-6 


(1.52 


± 0.34 


± 0.44 ) 


10" 


-5 


(1.33 ± 0.22 


± 0.20 )-10-^ 


n 


(MS) at 35.0 GeV 








(MS) at 38.3 GeV 






(MS) at 43.8 GeV 



1 (2.555 ± 0.009 ± 0.018)-10-^ 

2 (7.174 ± 0.052 ± 0.081)-10"^ 

3 (2.209 ± 0.026 ± 0.035)-10-2 

4 (7.42 ± 0.12 ± 0.17 )-10"* 

5 (2.694 ± 0.060 ± 0.091)-10-^ 



(2.509 ± 0.027 ± 0.018) 

(7.00 ± 0.15 ± 0.11 ) 
(2.154 ± 0.074 ± 0.064) 
(7.22 ± 0.35 ± 0.39 ) 
(2.61 ± 0.16 ± 0.22 ) 



10" 


1 


10" 


2 


10" 


2 


10" 


3 


10" 


3 



(2.371 ± 0.017 ± 0.026)-10-^ 
(6.316 ± 0.097 ± 0.098)-10-^ 
(1.884 ± 0.047 ± 0.034)-10-2 
(6.24 ± 0.22 ± 0.14 )-10-^ 
(2.26 ± 0.11 ± 0.07 )-10-^ 
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Table 4. Measurements of asiM^o) from event shape moments over the full analysed range of PETRA cm. energies, 14-44 GeV. 
The hadronisation uncertainty is taken to be the larger of the deviations observed using HERWIG and ARIADNE 





((l-T)i) 




(Bt> 






as(Mzo) 


0.1276 


0.1241 


0.1157 


0.1308 


0.1346 


Statistical error 


0.0004 


0.0003 


0.0002 


0.0004 


0.0009 


Experimental syst. 


0.0013 


0.0010 


0.0006 


0.0014 


0.0016 


HERWIG hadr. corr. 


-0.0017 


-0.0017 


-0.0003 


-0.0007 


+0.0011 


ARIADNE hadr. corr. 


+0.0002 


+0.0000 


+0.0009 


-0.0042 


-0.0051 


Hadronisation error 


0.0017 


0.0017 


0.0009 


0.0042 


0.0051 


Xf, variation: 












x^ = 2.0 


+0.0084 


+0.0076 


+0.0055 


+0.0097 


+0.0079 


x^ = 0.5 


-0.0068 


-0.0061 


-0.0043 


-0.0005 


-0.0059 


xVd.o.f. 


14.9/5 


16.7/5 


48.8/5 


98.8/5 


40.0/5 





{(l-Tf) 










(Mi) 


as(Mzo) 


0.1447 


0.1417 


0.1333 


0.1327 


0.1369 


0.1294 


Statistical error 


0.0008 


0.0005 


0.0004 


0.0006 


0.0019 


0.0004 


Experimental syst. 


0.0019 


0.0017 


0.0011 


0.0021 


0.0016 


0.0011 


HERWIG hadr. corr. 


+0.0009 


-0.0001 


+0.0006 


-0.0006 


+0.0026 


+0.0051 


ARIADNE hadr. corr. 


+0.0009 


+0.0007 


+0.0011 


-0.0048 


-0.0043 


-0.0024 


Hadronisation error 


0.0009 


0.0007 


0.0011 


0.0048 


0.0043 


0.0051 


Xft variation: 














x^ = 2.0 


+0.0141 


+0.0134 


+0.0125 


+0.0074 


+0.0088 


+0.0062 


x^ = 0.5 


-0.0113 


-0.0109 


-0.0103 


-0.0055 


-0.0067 


-0.0043 


xVd.o.f. 


13.5/5 


16.3/5 


33.7/5 


64.7/5 


13.7/5 


92.7/5 





{(1-T)'> 












as(Mzo) 


0.1514 


0.1497 


0.1434 


0.1376 


0.1352 


0.1364 


Statistical error 


0.0013 


0.0007 


0.0007 


0.0011 


0.0030 


0.0007 


Experimental syst. 


0.0022 


0.0021 


0.0014 


0.0032 


0.0027 


0.0012 


HERWIG hadr. corr. 


+0.0033 


+0.0016 


+0.0018 


-0.0006 


+0.0033 


+0.0069 


ARIADNE hadr. corr. 


+0.0016 


+0.0015 


+0.0012 


-0.0059 


-0.0039 


-0.0030 


Hadronisation error 


0.0033 


0.0016 


0.0018 


0.0059 


0.0039 


0.0069 


Xij variation: 














= 2.0 


+0.0166 


+0.0164 


+0.0162 


+0.0084 


+0.0084 


+0.0087 


= 0.5 


-0.0132 


-0.0131 


-0.0130 


-0.0063 


-0.0064 


-0.0067 


xVd.o.f. 


12.1/5 


16.5/5 


23.8/5 


43.8/5 


6.0/5 


66.9/5 





{(l-Tf) 












as(Mzo) 


0.1553 


0.1546 


0.1489 


0.1392 


0.1333 


0.1399 


Statistical error 


0.0018 


0.0009 


0.0010 


0.0019 


0.0045 


0.0010 


Experimental syst. 


0.0024 


0.0024 


0.0017 


0.0042 


0.0045 


0.0013 


HERWIG hadr. corr. 


+0.0051 


+0.0031 


+0.0030 


-0.0009 


+0.0034 


+0.0083 


ARIADNE hadr. corr. 


+0.0022 


+0.0022 


+0.0013 


-0.0068 


-0.0039 


-0.0036 


Hadronisation error 


0.0051 


0.0031 


0.0030 


0.0068 


0.0039 


0.0083 


X/j, variation: 














x^ = 2.0 


+0.0183 


+0.0185 


+0.0187 


+0.0083 


+0.0079 


+0.0094 


= 0.5 


-0.0145 


-0.0146 


-0.0148 


-0.0060 


-0.0060 


-0.0073 


xVd.o.f. 


10.9/5 


17.3/5 


17.3/5 


24.4/5 


3.2/5 


47.0/5 





({i-m 




{B't) 




((y2°3)=> 


(MS) 


as(Mzo) 


0.1580 


0.1586 


0.1525 


0.1397 


0.1314 


0.1411 


Statistical error 


0.0027 


0.0011 


0.0015 


0.0035 


0.0070 


0.0013 


Experimental syst. 


0.0029 


0.0025 


0.0020 


0.0052 


0.0061 


0.0017 


HERWIG hadr. corr. 


+0.0066 


+0.0044 


+0.0040 


-0.0013 


+0.0031 


+0.0094 


ARIADNE hadr. corr. 


+0.0027 


+0.0029 


+0.0012 


-0.0077 


-0.0043 


-0.0040 


Hadronisation error 


0.0066 


0.0044 


0.0040 


0.0077 


0.0043 


0.0094 


X/j variation: 














x^ = 2.0 


+0.0198 


+0.0204 


+0.0206 


+0.0078 


+0.0075 


+0.0096 


X/j, = 0.5 


-0.0155 


-0.0159 


-0.0161 


-0.0055 


-0.0057 


-0.0073 


xVd.o.f. 


9.6/5 


18.4/5 


11.9/5 


10.5/5 


17.3/5 


32.4/5 
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Table 5. Measurements of 03(^^20) from event shape moments over the full analysed range of PETRA cm. energies, 14- 
44 GeV, and the full range of LEP cm. energies, 91-209 GeV. The hadronisation uncertainty is taken to be the larger of the 
deviations observed using HERWIG and ARIADNE. The experimental systematic errors are estimated by the minimum overlap 
assumption 













as(Mzo) 


0.1271 


0.1242 


0.1165 


0.1259 


Statistical error 


0.0002 


0.0002 


0.0001 


0.0005 


Experimental syst. 


0.0008 


0.0006 


0.0005 


0.0016 


HERWIG hadr. corr. 


-0.0017 


-0.0020 


-0.0014 


+0.0006 


ARIADNE hadr. corr. 


+0.0023 


+0.0021 


+0.0019 


-0.0013 


Hadronisation error 


0.0023 


0.0021 


0.0019 


0.0013 


Xfi variation: 










Xu = 2.0 


+0.0078 


+0.0071 


+0.0053 


+0.0060 


Xu = 0.5 


-0.0063 


-0.0057 


-0.0041 


-0.0045 


xVd.cf. 


31.2/9 


36.1/9 


114/9 


173/9 





(il-Tf) 






(Bw> 


((2/?3)^> 




as(Mzo) 


0.1437 


0.1414 


0.1338 


0.1271 


0.1279 


0.1253 


Statistical error 


0.0005 


0.0003 


0.0003 


0.0004 


0.0012 


0.0003 


Experimental syst. 


0.0013 


0.0011 


0.0008 


0.0013 


0.0022 


0.0009 


HERWIG hadr. corr. 


+0.0007 


-0.0002 


-0.0002 


-0.0009 


+0.0018 


+0.0034 


ARIADNE hadr. corr. 


+0.0026 


+0.0026 


+0.0019 


-0.0021 


-0.0009 


+0.0003 


Hadronisation error 


0.0026 


0.0026 


0.0019 


0.0021 


0.0018 


0.0034 


Xu variation: 














Xu = 2.0 


+0.0131 


+0.0126 


+0.0122 


+0.0061 


+0.0066 


+0.0052 


Xu = 0.5 


-0.0107 


-0.0103 


-0.0100 


-0.0045 


-0.0051 


-0.0036 


xVd.o.f. 


22.3/9 


24.7/9 


47.0/9 


230/9 


46.4/9 


247/9 





{(l-Tf) 




(Bt> 






(Mi) 


as(Mzo) 


0.1509 


0.1495 


0.1436 


0.1317 


0.1282 


0.1307 


Statistical error 


0.0010 


0.0005 


0.0005 


0.0009 


0.0021 


0.0004 


Experimental syst. 


0.0016 


0.0014 


0.0011 


0.0020 


0.0025 


0.0012 


HERWIG hadr. corr. 


+0.0026 


+0.0011 


+0.0011 


-0.0008 


+0.0022 


+0.0047 


ARIADNE hadr. corr. 


+0.0027 


+0.0031 


+0.0018 


-0.0034 


-0.0012 


-0.0002 


Hadronisation error 


0.0027 


0.0031 


0.0018 


0.0034 


0.0022 


0.0047 


Xu variation: 














Xu = 2.0 


+0.0158 


+0.0154 


+0.0159 


+0.0069 


+0.0067 


+0.0071 


Xu = 0.5 


-0.0127 


-0.0125 


-0.0128 


-0.0052 


-0.0051 


-0.0055 


xVd.o.f. 


15.9/9 


22.0/9 


28.8/9 


117/9 


16.0/9 


194/9 





{(l-T)^) 












as(Mzo) 


0.1555 


0.1551 


0.1490 


0.1340 


0.1284 


0.1329 


Statistical error 


0.0015 


0.0006 


0.0009 


0.0017 


0.0036 


0.0007 


Experimental syst. 


0.0020 


0.0016 


0.0014 


0.0029 


0.0040 


0.0015 


HERWIG hadr. corr. 


+0.0043 


+0.0022 


+0.0023 


-0.0009 


+0.0024 


+0.0056 


ARIADNE hadr. corr. 


+0.0030 


+0.0037 


+0.0016 


-0.0046 


-0.0019 


-0.0008 


Hadronisation error 


0.0043 


0.0037 


0.0023 


0.0046 


0.0024 


0.0056 


Xu variation: 














Xu = 2.0 


+0.0179 


+0.0177 


+0.0185 


+0.0070 


+0.0067 


+0.0076 


Xu = 0.5 


-0.0142 


-0.0141 


-0.0146 


-0.0051 


-0.0051 


-0.0058 


x'/d.o.f. 


13.0/9 


21.7/9 


20.3/9 


50.2/9 


6.6/9 


139/9 





{(1-Tf) 










(M|> 


as(Mzo) 


0.1588 


0.1598 


0.1528 


0.1349 


0.1277 


0.1336 


Statistical error 


0.0024 


0.0007 


0.0014 


0.0032 


0.0058 


0.0010 


Experimental syst. 


0.0029 


0.0016 


0.0019 


0.0044 


0.0065 


0.0021 


HERWIG hadr. corr. 


+0.0059 


+0.0031 


+0.0034 


-0.0012 


+0.0024 


+0.0062 


ARIADNE hadr. corr. 


+0.0032 


+0.0042 


+0.0014 


-0.0056 


-0.0025 


-0.0013 


Hadronisation error 


0.0059 


0.0042 


0.0034 


0.0056 


0.0025 


0.0062 


Xu variation: 














Xu = 2.0 


+0.0197 


+0.0198 


+0.0206 


+0.0067 


+0.0065 


+0.0076 


Xu = 0.5 


-0.0155 


-0.0156 


-0.0160 


-0.0046 


-0.0049 


-0.0058 


xVd.o.f. 


11.6/9 


23.6/9 


13.9/9 


19.0/9 


3.1/9 


93.3/9 



